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Abstract 

Electrostatic bender optics are derived up to second order (third 
order in fields and the Hamiltonian) and applied to the proposed EDM 
proton ring. The results for linear optics agree with those already 
presented by V. Lebedev (Nov. 18, 2013). Second order optics is not 
sensitive to the shape of the fringe fields and formulas are given. It is 
shown that the proposed electrode shape that linearizes the vertical 
electric field is no advantage to this order. 


1 Introduction 

The general electrostatic bender has a given curvature in the median plane 
(arbitrarily assumed here to be horizontal), and, to control the motion in the 
vertical direction, also has curvature in that direction. It is assumed that 
the two curvatures are constant radius, but different, so the equipotentials 
are sections of a torus. Lebedev[l] has mentioned that the electrodes can be 
shaped in such a way that the vertical held is linear, i.e., and so has no xy 2 
coupling term. This shape will also be considered. 

Though the electrodes would be continuous in the vertical plane, they are 
necessarily truncated in the horizontal plane to allow the particles to enter 
and exit the held. Thus the radius of curvature along the reference trajectory 
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changes continuously from infinity outside the held region to a finite value in¬ 
side. The electric potential of such a truncated torus is derived to third order 
in deviations from the reference trajectory; the third order term contains the 
second derivative of the curvature function thus enabling analysis of lowest 
order fringe held effects. The relativistically correct Hamiltonian is derived 
also to third order. A canonical transformation is found which simplifies 
the Hamiltonian by eliminating derivatives of the curvature function. This 
allows an analysis of the linear optics and simple formulas for second order 
aberrations arising both from the curved geometry and the fringe helds. The 
linear optics agree with known; the second order aberration formulas were 
presented previously for the non-relativistic limit [2], 


2 The Potential 


In a general 
given by 


curvilinear system the Laplacian of a potential V(q l ,q 2 ,q 3 ) is 
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where the hi 


are the Lame coefficients and in general are functions of the q 1 .. 


2.1 Toroidal coordinates 

We start with a curvilinear coordinate system that conforms to the shape 
of the equipotentials: the coordinate along the intended reference orbit we 
call q 3 = s, the radially outward coordinate is q 1 = x, and the vertical is 
q 2 = y. The horizontal curvature (reciprocal of radius) is h(s), and the 
vertical curvature is k, a constant. In such a system, the Lame coefficients 
are hi = 1, h 2 = 1 + kx , h 3 = 1 + h(s)x. 

But in this system, the potential has no dependence on y, the direction 
orthogonal to x and s. We can therefore expand to third order as 1 

2 S 

__ rjQ ** 

H(a:,s) = oi 0 (s)x + a 20 (<s)— + a 30 (s)— (2) 

2 o 

1 This potential is dimensionless. Multiply by /3c(Bp) = ficpo/q = f3 2 'ymc 2 /q to get the 
potential. 
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We know the electric field on the reference orbit imposes the curvature h, 
so aio(s) = h(s). We can find the other coefficients from Laplace’s equation 
and eqn. 1: 


V(x, s ) = hx — \h(k + h)x 2 + \ 2 h (k 2 + kh + h 2 ) — h'‘ 
2 6 L v ' 


x 


(3) 


2.2 The Potential in Frenet-Serret Coordinates 

The Frenet-Serret coordinate system is the usual one of accelerator physics. 
It has no curvature in the ^/-direction, so the factor h 2 = 1 in eqn. 1. The 
potential in the median plane is the same in both systems, we need only 
add the following two more terms allowed by the symmetry: ao 2 (s)y 2 /2, 
a r2 (s)xy 2 /2. 

V(x, y,s) = V(x, s ) + i a 02 (s)y 2 + ^ a 12 (s)xy 2 . (4) 

We again use the curvilinear Laplacian (1), but this time with h 2 = 1, 
set it to zero and solve term-wise in the expansion to find a.Q 2 and a \ 2 : 
h Ji Ic Jc Ji 1 

V(x,y,s ) = hx— — (k+h)x 2 +— y 2 — — (2k+h)xy 2 +-[2h(k 2 +hk+h 2 )—h"]x 3 
2 2 2 6 

(5) 

Lebedev’s[l] linear vertical electric field (VL) potential also can be aug¬ 
mented to include lowest order fringe field effects and still satisfy Laplace to 
third order. The result is slightly different from the one above: 

Jl Jl /f 1 

V(x, y, s) = hx - -(k + h)x 2 + — y 2 + -[h 2 (k + 2 h) - ti']x 3 (6) 
2 2 o 

The reason for the difference is that potential 5 is derived for electrodes 
whose vertical sections are exact circular arcs while the potential 6 was de¬ 
rived with the specific intention to zero the xy 2 term. We shall see how the 
dynamics of the two compare. 


3 Dynamics 

3.1 General Hamiltonian for Electrostatics 

First, we must decide on the dynamical coordinates, especially the longitudi¬ 
nal (third) generalized momentum, which has been the source of considerable 


3 



confusion. A common approach, used especially by those accustomed to only 
magnetic elements, is to let the third momentum be A p/p. This is not the 
simplest approach, since when electric fields are included, it is not conserved. 
This means that when a particle enters the electrostatic element off-axis, it 
must receive a “kick” to get into the potential field. This kicks p, but leaves E 
unchanged. For an electrostatic bend of radius A, this kick is A p/p = Axj A\ 
the upper sign is for entry, the lower for exit. 

A better approach is therefore to let the third momentum be E. Since 
the fields are static, E is conserved; no kicks are required. 

For independent variable s, on a reference trajectory curving in the xs- 
plane with curvature h(s) = 1 /R(s), the Hamiltonian H = — p s is well- 
known, and 1 will not derive it here: 

H = -(1 + My( ^~ g$ ) -mV-pg-pg (7) 

We write E = E 0 + AE , and note that the large quantity under the square 
root sign is 


pi = Eq/c 2 — m 2 c 2 = (y 2 — l)m 2 c 2 = (/ 3ymc ) 2 , 
the square of the reference momentum. 


( 8 ) 


H 


■(1 + hx)p 0 


\ 


2Eq(AE — q<&) f AE — 

1 2 9 I 

Poc- V Po c , 


Po Po 


(9) 


Let us transform so that the third coordinate is not time, but a relative 
distance deviation w.r.t. the reference particle: i.e. from (t,—AE) to (r, p T ) 
where r = s — (3ct , p T = AE/(j3c). The generating function is 


F(t,p r ) = (A - f3ct)p r 

then the new Hamiltonian is K = H + dF/ds = H + p r 


( 10 ) 


K = Pr 


(1 + hx)p 0 


N 


1 + 2 



pepo) \Po 



Pi _p\ 
Pi Po 
( 11 ) 
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This cleans up considerably if we scale all momenta and the Hamiltonian 
by Po■ Px = Px/po, Py = Py/po, P-r = Pr/po, K = K/p 0 , and introduce the 
scaled potential V = 

K = P T -(1 + hx)yjl + 2 (P T -V)+ (3 2 (P T - V) 2 - PI - P 2 (12) 

This also makes the momenta accord with the more usual definitions, since, 
Px = x', P y = y\ and outside the electric field, P T = A p/p. 

This is the Hamiltonian of the electrostatic bench It is exact. 


3.2 Linear Optics 

Linear optics arise from second order Hamiltonian terms. Expanding to 
second order we get for both the potentials eqns. 5, 6: 

2 2 

V = hx - h(h + k)^- + hk V —, (13) 

where h = 1/Rq and k = 1 /R y is the curvature in the non-bend direction. 
The first term, needed to ensure that the reference trajectory is x = 0, yields 
the required electric field on the reference orbit: 


d$ _ _(3cpo dV_ _ {3^Eo 
dx q dx x=Q R 0 q ' 


(14) 


In the non-relativistic limit, the electric field is twice the beam kinetic energy 
divided by charge and bend radius: qS = mv 2 /R 0 . 

The first order terms in the resulting Hamiltonian all cancel as they should 
or the reference orbit would not be an orbit, so when expanded to second 
order it is 


p 2 P 2 

Ki = Pl + ^L + 
2 2 


El 

2y 2 


2 -E 

Rn 


xP T + 


2i?n 


X 2 + 


rp 

2R, 


2V 


(15) 


The parameters p and q are introduced as they parametrize the x and y 
focusing strengths: £ 2 + q 2 = 2 — (3 2 = 1 + y -2 , q 2 = k/h — R 0 /R y , R y being 
the curvature radius in the non-bend direction. (■ q 2 is given the symbol m in 
proton EDM group convention.) In the non-relativistic limit, for a cylindrical 
bend, £ = y/2, q = 0; for a spherical bend, £ = rj = 1. 
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The transfer matrix is easily derived from this Hamiltonian K\\ 


cos ^9 
Ro 


— jk sin £0 


0 

0 


2 -I3 2 


sin £9 


2 -13' 




^ sin £0 
cos £0 
0 
0 

-i?o(l — cos£ 0 ) 
0 


0 

0 

cos r rj9 

-ik sinr i d 

0 

0 


0 

0 


0 

0 


^ sin rtO 0 

77 ' 


cos )]9 
0 
0 


0 

1 

0 


^jf-R 0 {l ~ cos £0) 

s i n & 

0 

0 


R n 9 


1 ( 2 -py 

! 11 sin?®! 

p 2 v e ; 

V (0 J 


( 16 ) 

This matrix agrees precisely with that given by Lebedev[l], except that the 
56 element explicitly includes the velocity dependent part I/ 7 2 , allowing to 

read directly the phase slip factor for a solid electrostatic ring: ^ — (y~p) ■ 
Moreover, for such a ring, the horizontal tune is Q x = £ = y/2 — f3 2 — m, 
vertical tune Q y = rj = \fra, dispersion = 2 ^ 2 _ m Ro- 

1 have used this matrix to find the first order optics of a ring with 14 
benders with rf = m = 0.199 and separated by 0.834 metre drifts: tunes, 
^-functions, dispersion and phase slip factor agree precisely with the results 
presented by Valeri Lebedev[l], 


3.3 Second Order Aberrations 

Expanding to third order, the Hamiltonian for the toroidal electrodes be¬ 
comes 


K-2 = Ad - 


A 


f p 2 ' 

p2 + p2 _|_ £r 
± x ~ ± y ' o 
^ , 


+xy 2 — k 2 h + 


— P T x 2 
kh 2 ' 


kh h 2 ' 

+ 2—r 


, 27 ^ 7 

'A; 2 /r kh 2 


- PtV‘ 


kh 

v 


+ x P„ + P„. +2- 
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+ a: 
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For the potential 6, it is only slightly different: 
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' 3 \ /r 3 / 3 \ /i"V , 

v + 7 2 ) 6 ( 7 2 ) 6 A 7) 


p / p 2N 

K 2 = K 1 -^ [I# + Pj + -L 

1 j 


Pt-x 2 


' kh h 2 \ n 2 kh / 2 ~ \ , 

77777 + 2— - P T y —— + x[P x + P + 2— h 


w (i^!) +I 3 


2y 2 7 Z 

k 2 h kh 2 
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2y 2 
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Thus even though the latter potential was designed to eliminate the xy 2 
term, it re-appears in the Hamiltonian. The reason can be traced to the 
factor 1 + hx in front of the square root; coupled with the focusing term 
oc khy 2 this results in a term oc kh 2 xy 2 which is in fact larger than the term 
that was eliminated. Thus it is not clear that special non-circular-section 
electrode shapes are warranted. 


3.4 Fringe Field Shape Effect 

The presence of the second derivative of the curvature in the Hamiltonian 
seems to suggest that the effect in this order of the fringe field can be min¬ 
imized by shaping it. For example, if h(s) is linear, its contribution is zero. 
However, this technique proves to be futile: there necessarily will be regions 
of large h" at either end of such a linear fringe field, and these will have the 
same integrated effect as a fringe field with no linear part. To be sure, third 
(often called octnpole) and higher orders do depend on the fringe field shape; 
in general, the smoother and more extended the fringe field, the better. 

In fact, the second (often called sextupole) order effects of the fringe field 
are insensitive to the fringe field shape. This follows from the fact that the 
second derivative can be transformed away with a canonical transformation 
to a new horizontal transverse coordinate. 

We transform (x, P x ) —> (. X , P\) by means of the following generating 
function: 


,x 


x 2 P 


x 


G(x, P x , s) = xPx + h'-r-h 

o 2 

1 1 
X = x — -a; 2 /?. and P x = P\ — xPxh + - x 2 h' 


(19) 

( 20 ) 


The new Hamiltonian is K(X, Px, y, P y , z, P z ; s) = K + 


p2 

K = ^X 
2 


5 2 

2 


p 2 

+ ^ + x^+X 2 


2 h 
7 2 


P 2 X^ + P t 


2y 2 

P 2 
y 

2 2 


kh h 2 h 2 \ 2 kh 

T + ¥ + 2vJ +r T 


P T X h + 


h_ 

sy 2 


+ ( 21 ) 
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2y 2 


Xy 2 -k 2 h + 


kh 2 ' 


+ X 3 -k 2 h--kh 2 - 


kh 2 h 3 h 3 ' 
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1(23) 
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As (x, P x ) and (. X , Px) differ from each other only in higher order, the linear 
optics is unchanged. The Hamiltonian can be written as a sum of linear part 
K\ ( 21 ), momentum-dependent nonlinear part K P ( 22 ), and spatial nonlinear 
part K q (23). 


4 Results for Kq 


With no derivatives of h appearing, we can now to a good approximation 
assume h = l/Ro, a constant, and k = m/Ro . Further, scale lengths to be 
in units of R 0 . Then 


Kq = X 


(771 


1) 


2 1 
■f-2 

'*V Z 


I +Xy 2 -m 2 + 


m 

2j* 


= 0.792X 3 +0.024AV 
(24) 

The numerical values are for the proton EDM ring (7 = 1.2481, m = 0.199). 

The same transformation applies to the (VL) potential that is designed 
to eliminate the xy 2 term. In that case only Kq is different. The result is: 


Kq = X 3 


5m 1 m Vxy (| + J2Lj= 0.746X3 + 0.163*y 


"I-2 _ 

7 


6 ' 7 2 2y 2 


(25) 

As stated, the non-toroidal electrodes that make the vertical held linear 
appear to be of no advantage. 
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